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Phase sensitivity of perturbative nonlinear interactions 
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Despite the current concentration on phase control in few-cycle pulses, it emerges that there exists 
a wide class of nonlinear optical interactions in which the carrier phase is essentially irrelevant, even 
for the shortest pulse profiles. Most parametric processes and most perturbative processes fall into 
this category, although others such as above threshold ionization (ATI) do not. In an envelope 
approach, the carrier oscillations are not part of the problem because they are removed at the 
outset. When they are reinstated at the end of the calculation, one is free to include arbitrary phase 
shifts - within certain constraints. In many cases the constraints are relatively weak, and it follows 
that a single envelope solution can be used with an infinite range of choices for the carrier phase. 



I. INTRODUCTION 

We have investigated the role that the carrier phase 
plays in the propagation and interaction of few-cycle 
pulses affected by perturbative nonlinear interactions, us- 
ing theoretical methods as well as envelope approaches 
and numerical solutions of Maxwell's equations. The pur- 
pose of this paper is to examine in detail the process and 
consequences of how perturbative nonlinear interactions 
become phase sensitive, a subject usually passed over with 
only brief comment, or in rather general terms. 

Here focus our attention on x^") nonlinear processes 
and the interaction of multiple field components. Some of 
these, such as self phase modulation (SPM) are obviously 
phase insensitive. Others, are phase sensitive, such as in 
the experiments of Jones et.al. Q who use interference 
with a super-continuum generated from part of the input 
pulse, or Morgner et.al. |^ who use frequency doubling. 
Whether or not a x*^"-' process is carrier phase insensi- 
tive depends both on the nonlinearity and the number 
of field components involved; it is not guaranteed merely 
because it involves optical pulses only a few cycles long. 
The reason interactions involving many-cycle pulses are 
regarded as phase insensitive has little to do with the role 
of the carrier phase, but is instead due to the indistin- 
guishability of the many similar field oscillations in such 
pulses. Phase sensitive x^^^ interactions involving degen- 
erate parametric soliton interactions have also been dis- 
cussed ^. Many non-perturbative processes sugh as 
above-threshold ionization (ATI) (see e.g. Apolonski et. 
al.Q]) are also phase sensitive, however these are outside 
the intended scope of this paper. 

We show why some perturbative nonlinear processes 
are phase insensitive, but others are phase sensitive by 
reference to how is the x^""* processes are modelled, and 
whether the various frequency and nonlinear polarization 
terms are or are not significant. By reference to the pro- 
cedure used when an envelope approach is adopted, we 
systematize the procedure for generating models, and can 
therefore clarify how phase (in)sensitivity occurs. The 
key feature of an envelope approach ||^ is that the car- 
rier oscillations are removed at the outset and the pres- 
ence of the carrier is not part of the solution. In the 



many cases where the carrier phase is irrelevant, arbi- 
trary phase shifts can be included when the carriers are 
reinstated at the end of the calculation. It follows that 
a single envelope solution can include an infinite range 
of possibilities for the carrier phase. It is important to 
stress that this argument remains valid even if the pulses 
are so short that attempting to find a solution based on 
an envelope approach is impractical. 

In order to illustrate the results discussed in this pa- 
per, we compare two types of simulation. Firstly, we do 
envelope-based simulations of pulse propagation, which 
naturally agree with the phase insensitivity predictions of 
the envelope-based theory. Since we wish to emphasize 
that there is nothing intrinsically many-cycle or slowly- 
varying about these predictions, we propagate our pulse 
envelopes using the generalized few-cycle envelope ap- 
proximation (GFEA) theory @, ||, |9). This allows us 
(for example) to accurately model the self steepening 
of the pulse envelope in self phase modulation (SPM). 
Secondly, we do pseudospectral space-domain (PSSD) 
simulations [|l0) that solve Maxwell's equations without 
extra approximation. By varying the position of the 
carrier-like oscillations of the PSSD initial conditions, 
we are able construct a "pseudo-envelope" at each fre- 
quency of interest - the maximum excursions of a set of 
PSSD simulations of varying phases. For a set of results 
that match a phase insensitive envelope calculation, the 
pseudo-envelope should mimic the true envelope. 



II. BACKGROUND 

When using an envelope theory (e.g. SVEA or GFEA) 
to solve a pulse propagation problem, we find the carrier 
exponential (along with its phase (f)) often vanishes from 
the description, leaving only the chosen wavevector k and 
frequency lo as propagation parameters. The question 
we now ask is: Under what conditions can we apply a 
carrier with a different (p to the final envelope, and still 
reconstruct a physically correct answer? If we can change 
the underlying carrier phase, our propagation is phase 
insensitive, and for each envelope solution we can use 
an infinite set of phase choices. We will have solved for 
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an infinite set of input electric field pulses with a single 
envelope propagation. 

To examine this, we need to consider the definition of 
a pulse envelope and carefully consider the role of the 
carrier phase. An optical pulse can be written in the 
form 

with a carrier g^ikz-^ut+tp) ^^^^ envelope A(t,z). A 
shift in the phase of this pulse could be effected by 
either adjusting the carrier (adding the shift 9 to its ex- 
isting reference phase (jj) ; or by multiplying the envelope 
A{t,z) by e'^. We always choose the latter option, so 
when the phase of the pulse evolves during its propaga- 
tion, it is the phase of the envelope A that changes whilst 
the carrier remains fixed. If, for example, the centre of 
the pulse (where kz = tot) evolves from a cosine-like form 
at z to a sine- like form at z' , the sole effect is on the en- 
velope, i.e. A{t, z) — > A{t, z') = iA{t, z). 

There is nothing intrinsically many (or few) cycle 
about our conclusions here. If we can solve for the en- 
velope propagation in principle^ then our predictions for 
phase sensitivity will still hold. Of course some of the 
contributions neglected by the SVEA and GFEA theo- 
ries do depend on the carrier phase, but in general, for 
most materials, GFEA propagation is remarkably robust 
even for pulses containing only a few carrier cycles. 

The x^"'' processes we consider here model a wide va- 
riety of nonlinear optical behaviour. For these, nonlinear 
polarization affecting the pulse propagation of multiple 
field components is a simple product of the sum of those 
components 



P{t) - x^"^ 



[E^^-w]" (2) 

[E (^j(i)e'^'=^'""^*+*^'^ +C.C.)]" , (3) 



where kj , ujj , and 4)j are chosen to have convenient values, 
depending on the details of the system being described. 
When expanded, this will be a sum of many terms, each 
of which comprises a product of envelope contributions 
and carrier contributions. The normal procedure is to 
take each term in sequence, combining all its carrier con- 
tributions into a single exponential. We then match the 
time- like behaviour of this to one of the individual carrier 
exponentials jl^ - so the envelope associated with that 
carrier propagates according to the nonlinear contribu- 
tion that best matches its frequency of oscillation. This 
frequency match will often be exact because of the way 
we have chosen our set of carrier frequencies ujj , but usu- 
ally wavevector or carrier phase mismatches will remain. 
In most common situations, this carrier phase mismatch 
can be set equal to the phase of the frequency matched 
carrier, whilst still leaving some freedom in the overall 
choice of carrier phases. 

Note that the term "phase mismatch" is commonly 
used in nonlinear optics to denote the phase evolution 



resulting from a wavevector mismatch, and does not re- 
fer to a carrier phase mismatch. To avoid the potential 
confusion generated by this widespread terminology, in 
this paper we will only refer to either carrier phase mis- 
matches or wavevector mismatches, as appropriate. 

As we will see, some nonlinear effects impose no rela- 
tionship between the phases (pj of the field carriers; some 
allow a degree of freedom; and others leave no freedom 
at all. Any freedom to choose the <f>j in a particular case 
means that that case is phase insensitive. One important 
restriction on the ability of x*^"-* interactions to be phase 
insensitive is that they need to be (effectively) instanta- 
neous, so that the carrier oscillations at ujjt do not get 
folded into the calculation. For most nonlinear materials 
this is not a very significant restriction, as the nonlinear- 
ities are very fast compared to the optical frequencies of 
interest, but in the case of (e.g.) slower semiconductor 
nonlinear materials this may remain relevant. 

How these x'""* interactions stay phase insensitive (or 
not) will become clearer after looking at the specific ex- 
amples which follow. 



III. SIMPLE CASES 

A simple and commonly studied case of a nonlinear 
interaction is self-phase modulation (SPM). The under- 
lying nonlinearity that gives rise to SPM also gives rise 
to other terms, such as third harmonic generation (THG) 
terms A^) and other cross terms (e.g. ~ AiA*'^), Al- 
though SPM is phase insensitive, these other terms are 
not necessarily phase insensitive. In many cases, how- 
ever, they will be small and/or not wavevector matched 
- and if we are interested specifically in SPM, then our 
setup will be designed to ensure they are negligible. 

The polarization contribution that affects the field en- 
velope A during propagation under SPM is the "polar- 
ization envelope" Bspm, which has the same carrier as 
the field envelope. It is 



Bs 



PM 



X 



^^'>A.A*.A. 



(4) 



which has no dependence on the carrier phase. We 
demonstrate this in figjl , where a set of PSSD solutions of 
Maxwell's equations |1C| for SPM with a range of phases 
match up exactly with the phase insensitive GFEA enve- 
lope. The use of a GFEAQ rather than an SVEA pulse 
propagation equation for the envelope allows us to ac- 
curately reproduce the strong self-steepening that occurs 
for few-cycle pulses. It also emphasizes that this phase 
insensitivity still persists in the few-cycle regime. 

Another simple case is n-th order harmonic generation. 
This needs to be modeled using two field components, 
with frequencies ujn — noji. The underlying nonlinear- 
ity gives us a total polarization P{t) = x^^^E!{t)"', from 
which the two important contribution Bi, Bn affect the 
propagation of the field envelopes Ai , An respectively. As 
for the SPM case, other polarization contributions might 
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Comparison FDTD vs GFEA (BK Fig2o) (20040407) 



Comporison FDTD vs GFEA (x ) (20040419) 




FIG. 1: SPM: We match a GFEA envelope to a set of PSSD 
solutions to Maxwell's equations, and achieve excellent agree- 
ment. It is also possible to use the envelope to generate field 
profiles that match the individual PSSD solutions. 
Heavy line: GFEA envelope solution. 
Light lines: PSSD fields, phase intervals of 7r/4. 



in some general case be significant, but a harmonic gen- 
eration experiment would be designed to ensure they are 
negligible. The polarization envelopes Bi and Bn that 
match to the fundamental Ai and n-th harmonic An are 

Bi = x^"^^n^r"'exp{i[A„z + 0„-7^,/)i]}, (5) 
Bn - x^"^^>xp{z[-A„z + n(/)i -0„]}. (6) 

Here A„ = kn—nki is the wavevector mismatch. These 
have no dependence on ^^'s if <5 = — rKjii is fixed. This 
means there is a phase sensitivity, but it only fixes the 
relative phase between the field components. This leaves 
one free choice of carrier phase, and so the interaction 
can still be described as phase insensitive. 

To demonstrate the phase properties for harmonic gen- 
eration we present some simulation results for the special 
case of second harmonic generation, using two field com- 
ponents with frequencies uj2 = 2lui ; hence A2 = k2 — 2ki 
and the carrier phase condition is 02 = 2(/)i. In fig. ^ 
we show results with the initial condition A2 = 0, auto- 
matically satisfying the carrier phase condition, but still 
allowing a free choice of the carrier phase Conse- 
quently, it shows how a set of PSSD solutions (with range 
of initial (j>i ) closely match the phase insensitive GFEA 
envelopes. So, by ensuring that 02 = 20i, a single en- 
velope calculation can give us answers for any choice 0i 
value. 

In these simulations we needed to separate fundamen- 
tal and second harmonic fields (and hence their pseudo- 
envelopes) from the PSSD solutions for the total field. To 
do this we transformed into the spectral domain, sepa- 
rated out the fundamental and second harmonic spectral 
peaks, and transformed each back independently. This 
process involved some judgment; and also filtered out 
the low frequency optical-rectification effects, which were 
not included in the GFEA simulations in any case. This 
means that the resulting pseudo-envelopes seen on the 
graphs are not perfect matches to the envelopes, since, 
for example, the field spectra contributions in the region 
between the carriers will not shared out in exactly the 




FIG. 2: We match GFEA envelopes to a set of PSSD solu- 
tions to Maxwell's equations with excellent agreement. 
Heavy lines: GFEA envelopes Ai and A2; 
Light jagged lines: PSSD pseudo-envelopes ; 
Oscillating lines: PSSD second harmonic field E2 for (pi — 
(light), and the fundamental field Ei (dotted). 
The second harmonic envelope curve does not exactly match 
the pseudo-envelope because of both phase sampling and 
spectral filtering. The fundamental envelope and pseudo- 
envelopes barely be distinguished on the scale of this graph. 




FIG. 3: We match the fundamental and second harmonic 
GFEA envelopes to the corresponding PSSD solution to 
Maxwell's equations, and see excellent agreement (20i = 
02 = 0). 

Heavy lines: GFEA envelopes Ai and A2; 
Light lines: PSSD field split into Ei and E2. 
NB: peak _Bi(0) = 1.00 and peak £2(0) = 0.20. 



same way. 

To contrast with the above results, we now demon- 
strate how phase sensitivity can manifest itself in simu- 
lations that violate the 20i — 02 condition. We include 
a finite second harmonic pulse of varying phase 02 in the 
initial conditions, but fix 0i =0. On fig. |^ we show the 
case where 02 = 0, which is the one phase that satisfies 
201 = 02, and as a result we see the GFEA envelope and 
PSSD fields in agreement. In contrast, on fig. |^ we show 
the second harmonic PSSD fields for a range of 02, and 
we see that no sensible envelope function could reproduce 
the set of curves - as expected since 20i 7^ 02. 

In a set of PSSD simulations where both 0i and 02 
were varied, a given GFEA simulation would reproduce 
the PSSD data along a line where 20i — 02. Even if a 
particular case is phase insensitive, it is not guaranteed 
to be of the type we might find useful. 
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time (fs) time (fs) 

FIG. 4: Two slices of the PSSD second harmonic field E2 for 
a selection of 02 values with 201 7^ 02- In the wings of the 
pulse (left hand, 70-74fs) we can see approximately envelope- 
like behavior because the phase-sensitive polarization terms 
are small; but as we move towards the pulse peaks (right 
hand, 90-94fs), clearly no reasonable envelope function could 
match all the data. 

NB: peak EiiO) = 1.00 and peak £2(0) = 0.20. 

IV. COMPLETE EXAMPLE: THREE-FIELD x'^' 
INTERACTION 
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polarization contributions that are commonly considered, 
then step-by-step add the rest back in to the model, de- 
scribing the consequences for the phase sensitivity. 



A. Standard (resonant) case 

The standard subset of polarization contributions are 
the ones exactly resonant with one of the three carriers. 
As usual, Ba is a polarization envelope with the car- 
rier exp [i {kaZ ~ ujat + <^a)], SO that with the wavevector 
mismatch A = kp — ks — ki, we have 

B, = B,M = ApAy^'^+^'l'''-^''^^'') (8) 
Bs = BsM = ApA*e'^'+'^'l''''^^^-'*'''> (9) 

Each gives us (the same) carrier phase constraint (A): 
4>p = 4>s+4>ii and so any single envelope simulation can be 
turned into a multiplicity of solutions for the fields, using 
the two unconstrained carrier phases to determine the 
third. Thus, for resonant-only polarization contributions, 
the interaction is doubly phase insensitive (i.e. we have 
two free choices of carrier phase) . 



Phase sensitivity of PNL interactions 



For a three-field x^^'' interaction, such as in an optical 
parametric amplifier (OPA), we get a number of polar- 
ization terms to consider. In contrast to the previous 
section, which looked at common nonlinear polarization 
terms in isolation, here we give a more comprehensive 
review of all the nonlinear polarization terms in a partic- 
ular case. As we choose to use the interaction for an OPA 
or OPO system, we follow the usual naming convention 
and specify three field components (pump, signal, idler) 
with frequencies Wp > > uji such that ujp = ojg + uji. 
The polarization terms from P{t) = E{t)^ that affect the 
propagation of the field envelopes Ap, A^, Ai are - 



B. Wideband case 

Of course in a wideband (or miraculously well wavevec- 
tor matched) situation, other polarization terms might be 
relevant. Assuming that lOs — uJi ^ tOi, and 2uji ~ lUs, we 
would associate the extra polarization contributions with 
the propagation of the envelope with the nearest carrier 
frequency, so 



Bi — Bi^jf 
Bs = Bs^R 



Bi^w 
B 



w 



Bi^R 
Bs,R 



+ ^^^*g*('=^'-2fei)z+*('^'s-20i^21) 

(12) 



ks}z+i{24>i-<f>s) 



P+ = 



\2 i2{kpZ—ujpt-\- 
P 



+As 



{kiZ—uJit~{-4>i) 

^t{kpZ — UJpt'\-cf)p)j^*^~t{ksZ — LUst'\-(ps) 

^t{kpZ—ujpt'\-(pp) ^—t{kiz-ujit^(pi) _|_ j^^i2{ks 
^i{ks z—uJsi-\-(ps) A ^t(kiZ—u!it-\-(pi) 
^i{ks z—uJsi^(ps) A * ^-—i{kiZ—iAJit-\-(pi) 
^'i2{kiZ—L0it-\-(pi) 



-A 

-Ar. 



t+cf, ) Both of these add the extra phase constraint (B): (ps = 
2(j)i. Substituting this into constraint (A) above therefore 



gives us the new constraint 



3(/>i; which leaves us 



with just one unconstrained phase from which we can 
generate multiple field solutions from a single envelope 
-"=*+Mulation. 

For example, in some recent carrier-envelope phase sta- 
bilisation experiments |12, |l3|, the authors advance an 



The zero and negative frequency polarization compo- 
nents are not included here - the negative frequency ones 
(P~) apply to the evolution of the conjugate parts of the 



envelopes A* 



A* A*- 



and the zero frequency ones {P ) 



correspond to the uninteresting (for us) steady state con- 
tributions. 

In the following we will start by considering the three 



argument based (only) on the resonant terms to explain 
"4$!y they can generate an idler field insensitive to the ini- 
tial phase of the pump pulse (from which the signal is also 
derived). Since they have near resonant signal and idler 
fields, the two "wideband" terms above in cqns. (^ [T^ ) 
will produce components near the pump (since 2uji ~ ujp) 
and near the zero-frequency (since LOg — oJi ^ 0) respec- 
tively. Consquently their model is very insensitive to 
these contributions to the nonlinearity, especially when 
the pump field is strong. However, if they had chosen a 
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scheme with ujs—uji LOi, then the Ws — Wi nonlinear com- 
ponent would disrupt their phase independent idler gen- 
eration - although fortunately the term would be weak, 
since it is dependent on the idler strength. More serious 
problems would have ensued if the scheme required gener- 
ation of a signal pulse from pump-idler interaction, since 
then a strong oji + Ui term would disrupt the hoped-for 
generation of a phase-independent signal pulse. 



C. Extreme case 

Finally, in some extreme cases we might need (or want) 
to include the remaining polarization terms. Because of 
their high frequency, they are best assigned to drive the 



Up 



field component 



B„ — 



Bp,x 



(13) 



j^^j^^^i[{kp+ki)z~{ujp+u:i)t+<j>p+<pi] 

_j_^^^^g«[(fcp-|-fcs)^-('^p+'^s)t-|-0p + 0a] 

_|_^2gj2(fcs2:-i^st+0s) _|_ ^2 ^i2(kpZ-ujpt+4)p) ^^4^ 

These would add constraints 0i = 0, 0s = 0, 0„ = 



2(j)s, and 



respectively, any of which would leave 



no adjustable carrier phases. This means that a single 
envelope simulation corresponds to only one single field 
result. Trying to avoid this by adding extra carriers (e.g. 
by inventing a,ujq — 2ujs) is generally counter productive, 
as these will generate extra polarization terms, each of 
which will arrive with a new collection of constraints. 



V. HIGH HARMONIC GENERATION 

High harmonic generation is the process whereby an in- 
tense laser field is use to ionize an atom, with the electron 
recoUsion with the atom generating a wide range of high 
harmonics off the laser pulse. It is clearly a phase sen- 
sitive process, because for a short (few cycle) laser pulse 
it is possible that only one of the carrier-like field oscilla- 
tions is sufficiently near the peak of the pulse to reach the 
ionization threshold. Changing the relative position of 
the carrier w.r.t. the pulse peak could then shift the tim- 
ing of the above-threshold part of the electric field profile, 
and even possibly shift it from a positive value electric 
field to a negative value. This will either alter the timing 
or direction of the ejected electron, and hence its recoUi- 
sion and the generated HHG signal. The typical descrip- 
tion of HHG focuses on the electron trajectory, following 
a tunnel-ionization event (caused by the distortion of the 
atomic potential of the atom or molecule. Such models 
bear no relationship to the perturbative x^^'' models al- 
ready discussed here, and so cannot be interpreted in the 
way we found useful above. 

However, since HHG is a multi photon process, it is in- 
structive to imagine a high-order nonlinear "perturbative 



HHG" process which is compatible with our approach. 
We emphasize that here we are not trying to make an 
accurate model of the HHG process, but to reproduce 
some of its features in a way that enables us to relate 
its phase sensitivity to that of standard x^"^ nonlinear 
processes. 

In our simple model, we consider just one high-order 
multi photon polarization term, and show how the phase 
constraints it generates remove any possibility of phase 
freedom. Considering an odd high-order process, we 
specify N = 2m +1 (i.e. TV odd), and write down the po- 
larization terms resulting from a single envelope-carrier 
combination at the fundamental frequency w - 



Pati = -lE{t) 

N 



N 



(15) 



-7 



.71=0 



where C„ are the binomial coefficients. Thus the n-th 
polarization term drives the electron at a frequency of 
{2n-N)iu- 



P2n- 



N 



If we compare the effect of the n-th term on the elec- 
tron to that of the n'-th term, we see that they beat with 
an exponential dependence like 



^z(2n — 2n')(A;^ — cjt+(/ 



(18) 



Even this first comparison shows a clear phase sensi- 
tivity; before we have considered all A'' of the terms, and 
before we consider the likelihood of contributions from 
interactions of a different order (e.g. N'). 

In the cases discussed in previous sections, we could 
regard some terms as negligible because they were small 
in magnitude or poorly wave- vector matched. Here all 
terms contain N powers of the envelope amplitude, thus 
none can be dismissed as negligible in size. All, however, 
have different frequencies. Our comparisons above shows 
us that as long as the electron can respond significantly to 
even just two terms, it's response will be phase sensitive. 



VI. CONCLUSIONS 

We have shown that nonlinear interactions can run 
the full range from completely carrier phase insensitive 
through to completely constrained, as demonstrated for 
both the simple cases in section |l| and the detailed ex- 



amination of a three-field x^^-* interaction in section IV 



The principles and procedure outlined in section [V 



easily be applied to higher order x^"-* nonlinearities with 
differing numbers of field components - all that is neces- 
sary is a careful enumeration of the field components and 
significant polarization terms, followed by the process of 
matching up the polarization terms to the evolution of 
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the field components. Lastly, in section |V| we saw how 
any high-order x''"^ process will, without carefully de- 
signed conditions, be be phase sensitive purely because 
of the large number of constraints generated by the po- 
larization expansion. 

Although strictly speaking all nonlinear x^"-* interac- 
tions allow no residual carrier phase insensitivity, such 
complete descriptions are rarely necessary for accurate 
physical models. This is because in many realistic cases 
we can safely ignore many of the nonlinear polariza- 
tion contributions and so recover some phase freedom, as 
demonstrated in this paper. This is because the neglected 
terms will be far off resonance and/or not wavevector 
matched enough to accumulate and so will play a negli- 
gible part in the propagation. 

The phase (in)sensitivity of the particular models pur- 
sued in this paper relate to experiments involving the 
phase-control of few cycle pulses (e.g. ^) in the fol- 
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lowing way: a scheme will be phase sensitive, and hence 
useful to the experimentalists if, when analyzed, there 
is no freedom to alter the carrier phase under some pre- 
dicted final-state pulse envelope. This contrasts with the 
case for processes like SPM, where a single envelope prop- 
agation code can predict a final-state envelope where any 
carrier phase can be used to generate a valid final-state 
electric field profile. 

The phase insensitivity we have discussed can manifest 
itself over the complete range from many- to few-cycle 
pulses - as long as in principle an envelope solution to 
the propagation is possible - and in principle, we might 
even retain all of the terms approximated away in the 
GFEA pulse propagation theory. Our conclusions are 
supported by successful comparisons between envelope 
predictions of carrier phase insensitivity and sensitivity 
with PSSD solutions of Maxwell's equations. 
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